Realistic neutron star models in f(T) gravity

We investigate the nonrotating neutron stars in $f(T)$ gravity with $f(T)=T+\alpha{T}^2$, where $T$ is torsion scalar teleparallel formalism of gravity. In particular, we utilize SLy and BSk family equations state for perfect fluid to describe stellar matter search effects modification on models stars. For positive $\alpha$, results a smaller mass comparison general relativity, while will contain larger amount negative $\alpha$. Moreover, there seems be an upper limit central density $\alpha>0$, beyond which effective would have steplike phase transition pressure profiles, collapsing numerical system. obtain mass-radius relations realistic subject them joint constraints from observed massive pulsars PSR J0030+0451, J0740+6620, J2215+5135, gravitational wave events GW170817 GW190814. star model able accommodate all mentioned data, parameter $\alpha$ needs than $-4.295$, $-6.476$, $-4.4$, $-2.12$ (in unit ${G}^2M_\odot^2/c^4$) SLy, BSk19, BSk20, BSk21 state, respectively. If one considers unknown compact object event GW190814 not hence excludes this dataset, can loosened $\alpha<-0.594$, $-3.5$, $0.4$ $1.9$ ${G}^2M_\odot^2/c^4$),